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ABSTRACT 


\ 


A theoretical  analysis  is  presented  for  treating  the 
free  vibrations  of  submerged,  ring-stiffened  cylindrical 
shells  with  simply  supported  ends.  The  effects  of  the  ec- 
centric stiffeners  are  averaged  over  the  thin-walled  iso- 
tropic cylindrical  shell.  The  energy  method  is  utilized  and 
the  frequency  equation  is  derived  from  Hamilton's  Principle. 
All  three  degrees  of  freedom  are  considered.  Numerical  re- 
sults are  presented  for  the  frequencies  of  several  example 
shells.  Comparisons  with  previous  theoretical  and  experi- 
mental results  indicate  reasonably  good  agreement  for  shells 
immersed  in  water.,  but  poorer  agreement  in  vacuum.  The  in- 
accuracies are  due  to  limitations  of  the  Donnell  type  ortho- 
tropic shell  theory  used.  The  procedure  appears  well  adapted 
to  preliminary  optimization  studies  of  shells  immersed  in 
water  under  minimum  natural  frequency  constraints  and  for 
optimal  separation  of  the  lower  frequencies. 
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CHAPTER  1 


INTRODUCTION 

Since  the  early  efforts  of  Rayleigh  [ 1]  1 and  Love  [2], 
the  vibration  of  shells  in  a vacuum  has  been  extensively 
analyzed.  Typical  of  recent  work  on  prestressed,  eccentric- 
ally stiffened,  cylindrical  shells  of  finite  length  are  the 
investigations  by  McElman,  et  al.  [3],  using  the  orthotropic 
or  smeared  stiffener  approach  and  Harari  and  Baron  [4],  who 
consider  discrete  stiffeners. 

The  dynamic  interaction  between  shells  and  fluids  has 
also  received  considerable  attention.  Junger  [5,  6,  7]  was 
the  first  to  analyze  the  free  and  forced  vibrations  of  a 
cylindrical  shell  submerged  in  an  acoustic  medium.  He  treated 
an  infinitely  long,  cylindrical  shell  utilizing  plain  strain 
analysis.  The  transient  response  of  a submerged,  infinitely 
long,  ring- stiff ened  cylindrical  shell  has  been  studied  by 
Herman  and  Klosner  [8]  and  Lyons,  et  al.  [9].  A submerged 
cylindrical  shell  of  infinite  length,  having  radial  surface 
motion  over  a stiffened,  finite  section  has  been  studied  by 
Paslay , et  al . [ 10 ] . 

Recently,  some  interest  has  been  directed  toward  optimi- 
zation of  submerged  shells  with  eccentric  ring  stiffeneres 
under  natural  frequency  constraints  and  or  with  dynamic 
merit  criteria  [11-13].  Unfortunately,  this  early  work 


lumbers  in  brackets  designate  references. 
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utilizes  an  invacuo  vibration  model  which  yields  excessively 
large  frequency  values. 

The  objective  of  this  ■ ffort  is  to  develop  a more  ac- 
curate model  for  such  optimal  design  studies.  Donnell  type 
orthotropic  theory  along  the  lines  of  Ref.  [3]  is  used  in 
this  preliminary  study  because  of  its  computational  simplicity. 
Conservation  of  computational  effort  is  quite  important  for 
preliminary  optimization  investigations,  since  it  is  desir- 
able to  perform  parametric  optimal  design  studies  so  as  to 
determine  the  characteristics  of  optimal  designs  and  since 
optimization  for  even  a single  set  of  design  parameters  often 
requires  several  hundred  to  several  thousand  sets  of  fre- 
quency determinations.  Thus,  it  is  desirable  to  use  the 
simplest  model  capable  of  generating  reasonable  results  for 
such  studies. 

This  report  therefore  describes  a procedure  for  the  nat- 
ural frequency  determination  of  submerged,  eccentrically  ring 
stiffened,  cylindrical  shells.  The  effects  of  eccentric 
stiffeners  are  averaged  over  the  thin-walled  isotropic 
cylindrical  shell.  The  case  where  axial,  tangential,  and 
radial  motions  of  the  structure  are  considered  as  well  as 
where  only  radial  motions  are  included.  Boundary  conditions 
are  established  and  satisfied  at  the  fluid-structure  inter- 
face, yielding  expressions  for  the  radiated  pressure.  The 
derivation  of  the  frequency  equation  is  accomplished  by  de- 
fining the  Donnell  type  of  nonlinear  strain  displacements 
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for  the  shell  and  its  stiffeners,  formulating  the  potential 
and  kinetic  energies  of  the  system,  and  then  applying  Hamil- 
ton's Principle.  A 3t  of  appropriate  nonlinear  equilibrium 
equations  and  boundary  conditions  are  derived.  The  nonlinear 
equations  are  used  to  derive  linear  equations  that  govern 
small  vibrations  of  the  system  and  a three-degree  freedom 
and  a one  degree  of  freedom  frequency  equation  are  obtained. 
The  frequencies  are  numerically  obtained.  For  the  three 
degrees  of  freedom  formulation  once  the  values  of  frequency 
have  been  obtained,  the  corresponding  ratios  of  radial,  axial, 
and  tangential  amplitudes  may  be  evaluated.  The  numerical 
results  obtained  from  the  present  analysis  for  several  cases 
of  interest  are  given  in  Chapter  5.  Reasonably  good  agree- 
ment is  found  with  other  theoretical  investigations  and 
with  experimental  results  for  the  case  of  shells  immersed  in 
water.  Poorer  agreement  is  found  for  the  case  of  shells  in 
vacuo  or  air. 


- 
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CHAPTER  2 

GOVERNING  FIELD  EQUATIONS 

The  well  known  field  equation  for  a homogeneous  acoustic 
fluid  medium  [14]  is 


Written  out  in  terms  of  the  notation  used  herein: 


ili  + _i_  H + a 2<i>  + iii  = _j_  lii  (2.1) 

3r2  r ar  r2  3(y/R)2  3x2  c2  at2 

in  which  q is  the  velocity  potential  function,  c is  the 

velocity  of  sound  in  the  fluid,  t is  the  time  variable,  x 
is  the  shell  coordinate  parallel  to  the  axis  of  the  struc- 
ture, (Figure  2.1),  r is  the  variable  describing  the  distance 
from  the  shell  surface  into  the  fluid,  and  y is  the  circum- 
ferential coordinate  at  the  mean  shell  radius  R. 

The  velocity  potential,  <p , and  the  radiated  pressure, 

P can  be  written  as 
r 


<)> 


iwt 

- e 


4>cos 


_J!X 

R 


s in 


mux 

L 


(2.2) 


Pr  = ela)t  Pr  cos-nf-  sin-21^-  (2.3) 

in  which  $ is  the  amplitude  of  the  velocity  potential,  Fr 
is  the  amplitude  of  the  radial  pressure,  w is  the  circular 
frequency,  m is  the  number  of  axial  half  waves,  R is  the 
mean  radius  of  the  cylindrical  shell,  and  L is  the  length 
of  the  cylindrical  shell. 
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Substituting  equation  (2.2)  into  the  governing  field 

equation  (2.1)  for  the  fluid  medium  yields: 

9 2 $ + 1 9$  + (-  u2  . n2  _ m2-rr2  ^ $ = Q 

3r2  r 3r  c2  r2  L2 


Let 

K 2 - ( -3-  )2'  ( )2  (2.4) 

If  the  fluid  is  assumed  to  be  infinitely  extended,  the 
solution  of  equation  (2.4)  for  an  outgoing  wave  is  given  by 

<t>  = AHn  2 (Kr)  (2.5) 

in  which  H is  the  Hankel  function  of  the  second  kind  of 

n 

order  n,  H *-2^  ~ J - iY  . 

n n n 

The  constant  A is  evaluated  by  ensuring  that  the  velocity 

of  the  shell  and  the  velocity  of  the  fluid  medium  are  equal  at 

the  shell-fluid  interface,  r = R,  i.e. 

3w  _ 3<f>  (2.6) 

3t  = ~ Sr- 

For  a shell  simply  supported  at  each  end,  the  radial 
deflection  is  assumed  to  be 

w = Welwt  sin  cos  (2.7) 


in  which  W = the  amplitude  of  the  radial  deflection. 

Substituting  equations  (2.2)  and  (2.7)  into  equation 
(2.6)  results  in 


iajW  = 


3 $ 
3r 


(2.8) 


Upon  substituting  equation  (2.5)  into  equation  (2.8), 

A = — (2.9) 


the  constant  A can  be  written  as 

iwW 


KH'  ^ (KR) 
n 
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- 
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CHAPTER  3 

ENERGY  OF  THE  STIFFENED  STRUCTURE 

Consider  an  axisymmetric  structure  consisting  of  a 
cylindrical  shell  with  reinforcing  ring  stiffeners.  In  the 
cylindrical  portion  of  the  structure,  the  displacements  of  the 
shell  are  defined  by  the  three  orthogonal  components  u,  v, 
w,  in  the  x,  y,  z directions  respectively,  which  are  functions 
of  the  coordinates  x and  y only  at  the  middle  surface  (see 
Figure  2.1). 

Potential  Energy 

The  strain  energy  of  the  unstiffened  thin-walled  iso- 
tropic cylinder  [15]  is 


h 2ttR  L 


di  o o 

2 


+ 2lJexTeyT  + YxyT2  )dxdYdz  C3-1) 

in  which  e^,  and  YXyj  are  the  total  normal  and  shearing 

strains,  E is  Young's  modulus,  and  y is  Poisson's  ratio.  A 
cylinder  may  be  considered  thin  walled  when  the  thickness  h 
is  sufficiently  small  compared  to  the  radius  R.  Usually, 

R/h  > 60. 
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The  Donnell  type  nonlinear  strain  displacement  relations  used 
here  are 


ExT 

UT,x 

+ i w 2 

I w ’ X 

eyT 

vT,y 

+ w/R  + y 

w,  2 

y 

Vt  “ 

,y 

+ vT,x  + W 

’x  w,y 

(3.2) 

where 

* u - zw,x 


vT  = v - zw,y  - (3.3) 

The  quantities  u,v,  and  w are  the  displacements  of  the 
middle  surface  of  the  cylinder  wall.  Thus, 


’xT 


u’x  * \ v’x 


zw, 


XX 


yT 


v,y  + 5 + I W4 


zw. 


yy  (3.4) 


^xyT  * u’y  * T-x  * w’x  ' 2zw-xy 


After  integTat-ion  with  respect- to  z-,  the  cylinder  strain 


energy  becomes 


♦ 2 (1-u)  w,2xy  ] dxdy 

where 

cx  ■ exT  U - 0 

*y  ' VT  I x ’ 0 

xy  ' Vt  Iz  * 0 C3.6) 

If  the  displacement  in  the  cylinder  and  stiffening  rings 
are  continuous  and  the  properties  of  the  stiffening  rings 
are  averaged  over  the  spacing  \,  the  total  strain  energy  [11] 
for  the  stiffening  rings  of  spacing  \ attached  to  the  shell 
is  found  to  be 

J 


(3.S) 
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The  first  term  in  equation  (3.7)  is  the  strain  energy 
of  extension  and  bending.  The  quantity  dAr  is  an  element  of 
cross-sectional  area  of  the  stiffening  ring,  and  Gr<Jr  is  the 
twisting  stiffness  of  the  ring  section.  Substitution  of  the 
first  of  the  equations  (3.4)  into  equation  (3.7)  and  integrat- 
ing over  the  area  of  the  stiffening  ring  yields  the  following 
expression  for  the  stiffening  ring  strain  energy 
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Here  z"r  is  the  distance  from  the  middle  surface  of  the 
isotropic  shell  to  the  centroid  of  ring  cross-section,  and 
I is  the  moment  of  inertia  of  the  stiffening  ring  with  re- 
spect to  an  axis  in  the  middle  surface  of  the  isotropic 
shell.  It  should  be  noted  that  z is  positive  for  a stif- 
fening ring  on  the  outer  surface  of  the  shell  and  negative 
for  a ring  on  the  inner  surface. 

The  potential  energy  of  external  pressure  and  an 
externally  applied  axial  load  resultant  (positive  in 
compression)  is 


Pwdxdy  + 


o 


L 

o 


dy 


(3.9) 


I 


where 


function  of  in  this  case 


Ph  * Pr 


constant  external  pressure 


radiated  pressure 


The  quantity  e is  the  distance  from  the  middle  surface  of 
the  isotropic  shell  to  the  line  on  which  the  load  resultant 


5Tx  acts 


The  potential  strain  energy  of  the  combined  structure 
can  be  written  as  the  sum  of  the  energies  of  the  cylindrical 
shell,  the  stiffening  rings,  and  the  loads  as  follows: 


Vc  + Vr  + VP 


(3.10) 


Kinetic  Eners 


The  kinetic  energy  of  the  system  can  be  written  in 
terms  of  the  kinetic  energies  of  the  -cyl-inxi-rical  shell 
segments  and  the  stiffening  rings. 
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The  kinetic  energy  of  the  unstiffened  thin-walled 


cylinder  is 


p h (u2  + v2  + w2)d*dy(  3. 11) 


in  which  Pc  is  the  density  of  cylinder  material. 
The  kinetic  energy  of  the  stiffening  ring  is 


— (u2  + ^2+  (3 . 12) 

l 


in  which  p r is  the  density  of  ring  material. 


The  total  kinetic  energy  of  the  system  is 


J7  (ii2  + tf2  + w2)d*dy  (3.13) 


.i . . _ i. 


16 


> 

I 


CHAPTER  4 

FREQUENCY  EQUATION  DERIVATION 

Equation  of  Motion  and  Boundary  Conditions 

The  partial  differential  equations  and  boundary  con 
ditions  are  derived  from  Hamilton's  principle 


[T  - V]dt  = 0 


(4.1) 


where  T and  V are  the  structural  system’s  kinetic  energy 
and  potential  energy  respectively. 

The  three  motions  of  this  conservative  system  from  a 
given  initial  configuration  to  a given  final  configuration 
in  a time  interval  (t^,  t2)  are  obtained  by  allowing  variation 
of  the  three  displacements  <5u,  <5v  and  <5w  to  be  arbitrary 
and  utilizing  the  fundamental  technique  of  the  calculus  of 
variations  [16].  From  equation  (4.1),  the  three  equations 
of  motion  are  derived  as  follows: 


k 


R 


+ 


2 


7 w’xJJw’ 


+ 11  u’x 


+ 


yy 


ErAr 


Eh 


2(l+d) 


V,  + V,  W,  W,  >1 

y X ’X  ’y  J 


W, 


xy 


Mw  = 0 


C4.2) 


The  necessary  and  sufficient  conditions  to  cause  the 
equation  (4.1)  to  vanish  during  the  variational  process  (3.1) 
are  to  zero  the  remainder  part  in  addition  to  the  equation  of 
motions.  Then  the  natural  boundary  conditions  are  established. 
These  conditions  at  each  end  of  the  stiffened  cylindrical  shell 
are : 


D [ 


+ 


(2-U)  w,  1 - 


Eh 


f u. 
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N +■  M 
X , X wxy , y 


N 

xy,x 


Rv 


0 


- M + M 

x , xx  xy , xy 


M 

yx.xy 


M 


y»7 y 


- N w - N w - 2N 
X ™ ’ xx  y ’yy 


w 


xy  ’xy 


+ P 


Rw 


(4.4) 


The  boundary  conditions  (4.3)  may  be  rewritten  in 
terms  of  stress  resultants  as 


M 


x , x ( Mxy , y " Myx , y)  + Nx  w ’ x + Nxy  w ’ y 


or  w = 0 


Mv  + N e =0 

X X 


°r  w>x  = 0 


M = -(■ 
■ yx  \ 


Gh3 


GJ 
r r 


-)  w. 


(4.6) 


Vibration  Equations 

The  equations  of  motion  (4.4)  derived  in  the  last  section 
are  used  to  obtain  linear  equations  which  govern  the  small 
amplitude  vibrations  of  a prestressed,  eccentrically  stiffened, 
cylindrical  shell  which  is  submerged  in  a fluid  medium. 


The  deformations  associated  with  the  vibration  of  a 
prestressed  cylinder  are  divided  into  two  parts  as  follows: 


U = UA  + UB 


V = \ + VB 


w =»  wA  - wB 


(4.7) 


The  first  part,  denoted  by  subscript  A,  is  an  axisymmetric 
static  prestress  deformation  which  occurs  prior  to  excitation' 
at  one  of  the  natural  frequencies.  The  second,  part,  denoted 


by  subscript  B,  is  a small  additional  deformation  which 
occurs  as  a result  of  the  excitation.  Since  the  A sub- 
scripted quantities  are  static  they  are  axisymmetric 
deformations;  therefore,  the  terms  Mii Mv^,  MwA  and  all 
derivatives  with  respect  to  y vanish.'  The  equilibrium 
equations  which  govern  these  deformations  are  found  from 
equations  (4-4)  as 


"xA.x  * » 


N 


xyA,x 


NyA 

•M  , + -i-u 

xA,xx  R 


NxAwA,xx  * -Ph  * 0 


(4.8) 


If  there  is  no  applied  shear,  equation  (4.8)  yields 


0 


If  we  let  all  derivatives  with  respect  to  y equal 
zero  for  axisymmetric  deformation,  a set  of  appropriate 
boundary  conditions  are  found  from  equations  (4.5)  to  be 
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The  equilibrium  equations  and  boundary  conditions 
which  govern  the  dynamic  deformations  (subscript  B)  are 
obtained  by  substituting  equation  (4.7)  into  equation 
(4.4)  and  (4.5).  By  eliminating  the  axisymmetric  pre- 
stress equations,  and  retaining  only  linear  terms, 
i.e.  neglecting  the  high  order  terms  under  the  small 
amplitude  vibration  assumption,  the  following  equations 
governing  the  dynamic  deformations  (subscript  B)  are 
obtained: 


NxB,x  + NxyB,y  * = 0 


NyB,y  + NxyB,x  ' Mv3 


0 


■D  (wB,yy  * “ wB,xx) 


Vr 


B,yy 


EA 


+ z — ( v„  + — z W \ 

r i \ B,y  R r B,yy) 


Vb  * HT~  WB,xy 


_ /Gh  . r r \ w_ 

^y  xB  " ( 6 i ) B , xy 


(4.13) 


Frequency  Equation  Derivation 


Assuming  a constant  prestress  deformation  W^,  the  solution 


equations  are 


N . * - IT 

xA  x 


V - - - - PR 


(1.14) 


Equation  (4.11)  may  now  be  written  as 
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NxB,x  * NxyB,y  ' M UB  ' 0 


V,y  * NxyB,x  ‘ M *B  * 0 


-Mr,  +M„  "Mr,  M „ 

xB,xx  xyB,xy  yxB,xy  yB,yy 


N. 


* -I5-  * Vs, XX  * Vs.yy  * P1 


JT  v?B  * . 0 


(4 


The  displacement  ufi,  Vg  and  wg  which  satisfy  simple 
support  boundary  conditions  are  given  as 

u . v .it  cos  (J-L.)  si,  (JSjJ.) 


».  - ? eiut  sin  (-ap-)  sin  (-SJ-j 


w 


B 


W eiaJt 


(^eH  cos(-^-) 


15) 


sm 


(4.16) 
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The  equation  C2.ll)  gives 

Pr  3 Pw  u)2f(<*0  elut  W sin  cos  (— -fc-) 

where  m is  the  number  of  axial  half-waves  and  n is  the 
number  of  circumferential  full  waves.  If  equations  (4.16) 

are  substituted  into  equation  (4. 15),  taking  into  account  (4.13), 
and  the  immediately  preceding  equation,  it  yields 

. (W  ♦ -Ik sill  * Gh  2) ) u . 

V 1-y2  L R ' 

.*  Sr  (“Pu  * Gh)  7 ' ““j.  Sr  w ' 0 


mn-rr  /UEh 
LR  v 1-u 


+ Gh)  U 


Eh nf_ 

1-U2  R2 


+ 


ErAr 


n2 

R2 


Gh 


m2m2  \ t? 


1 TT  \ 

L2  ' 


n / Eh 
R2  1-U2 


+ 


E A 
r r 


+ 


ErVrn: 

n. 


o 
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uEh  umir  u n / 

1-U2  R2  ' 

E A 

Eh  r r 

1-U2  1 

E A z n 

+ r r r >17  - 

“ 

1R  / 

^ /GhJ  GrJr  \ mmr  2 ^ ErXr  . n 4 

x ) [tx  ) i _ (it) 

* ^ f1  * r 

1 R2  r 

/ 2n  2 A — n4  o Eh  1 

(-*-*'*  R*  R2 

1 

hi 

l 

|)  - [ST  + pwf(a))]tu2  17  = 

C 4. 17) 


The  last  of  these  equations  contains  a term  involving 
f(w)  which  may  be  real  or  complex,  depending  on  whether  the 
K2  of  Equation  (2.4)  is  negative  or  positive,  respectively. 
When  f Cu)  is  complex,  the  real  part  may  be  considered  as  the 
"added  mass"  effect  of  the  fluid  and  the  imaginary  part  a 
form  of  damping  associated  with  energy  loss  to  the  fluid. 

For  the  range  of  frequencies  of  interest  here  for  shells 
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immersed  in  water,  K2  , will  usually  be  negative,  and  thus 
f(to)  is  real.  Thus,  for  the  purpose  of  this  study,  a fre- 
quency equation  where  f(co)  is  real  is  developed. 

Equation  (4.17)  may  conveniently  be  written  as 

[A]  {B}  = 0 (4.18) 


where 


(B) 


(4.19) 


This  is  not,  however,  a linear  eigenvalue  problem,  since  to 
occurs  in  f (to)  . 

For  the  case  where  f(io)  is  real,  the  determinant  of 
the  displacement  coefficients,  when  set  equal  to  zero,  yields 
a sixth-order  equation  of  the  circular  frequency 

£a]  = 0 (4.20) 

where 

a _ , _ mnir  r ru  ± uE h 

Af  ->  = A^  - — (Gh  + ) 

1 -y 2 


i ji\ 


• ' *■) 


The  result  of  [A]  = 0 is 

A^ui  6 + + A^u)2  + Aq  = 0 (4.Z2  ) 


where 


r. 
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f7  ■ D [(T-)'  * (l>  * 2“(rF)2] 


and 


/ Gh'  r r \ ^ r i / “ N 

Ctt  “i — ; (nr)  i (r) 


3 GrJr  \ / mnir  2 


Vr  , n ‘ 


E_A. 


* ^ r * *r  *>] 


Eh  1 


1-u2  R2 


(=H 


K (IV 


(4.  23) 


FjFI2 


A2  - B(f1(F2  * F.)  - StFy] 


•!  - Fl[F2  f3  - CF4)!|-  H [(f2  +f3)  f7  - CF5)!-CF6)2] 
>0  =-[(F4)2  - F2  Fj]F7  * CFS)2  F2  * CFS)2F3 


+ 2F4  F5  F6 


(4.  24) 


..  . . 
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{ 

After  the  natural  frequency,  w,  is  obtained  by  equating 

is 

the  determinant  [A]  to  zero  and  substitution  of  the  the  value 
of  w into  equation  (4.18),  the  non-trivial  solution  will  pro- 
vide the  amplitude  ratios  from  the  three  algebraic  equations 
of  the  coefficients  of  U,  V AND  17. 

Where  only  one  degree  of  freedom  is  to  be  considered, 
only  the  radial  surface  motion  of  the  system  is  permitted 
and  other  surface  motions  are  neglected.  The  frequency 
equation  is  derived  from  equation  (4.20)  by  equating  Rco2 
terms  to  zero  in  A.^  and  f<22'  This  equation  has  the  fol- 
lowing form: 

A'  w2  + A"  = 0 (4.25) 

Numerical  Procedures 

Both  the  three  and  one  degree  of  freedom  frequency  equa- 
tions (4.22)  and  (4.25),  respectively,  must  be  solved  numer- 
ically since  some  coefficients  are  functions  of  co . For  this 
report,  a sequence  of  values  of  w were  substituted  into  the 
left  hand  side  of  equation  (4.22)  and  a bisection  method 
employing  the  second  rule  [17]  was  used  to  locate  a root  when 
a change  in  the  sign  of  the  left  hand  side  occurred.  This 
technique  was  used  to  investigate  the  possibility  of  multiple 
real  roots.  None  were  found.  A similar  procedure  was  used 
for  the  single  degree  of  freedom  equation  for  programming 
convenience . 


Such  procedures  are,  of  course,  too  computationally 


wasteful  for  optimization  work.  In  the  optimal  design  situa- 
tion, one  almost  always  has  a very  close,  or  may  obtain  a 
reasonably  close,  estimate  of  the  value  of  co  satisfying  equa- 
tion (4.18)  or  (4.25).  Thus,  one  may  effectively  use  suc- 
cessive substitutions  to  solve  these  equations  in  relatively 
few  iterations.  For  example,  equation  (4.18)  may  be  treated 
as  an  eigenvalue  problem,  where  the  f(w)  term  is  computed 
using  the  close  approximation,  to  obtain  a still  closer  ap- 
proximation. This  procedure  would  be  repeated  until  conver- 
gence occurred.  Equation  (4.25)  could  be  solved  similarly. 
Thus,  although  the  computational  effort  associated  with  fre- 
quency determination  of  submerged  structures  would  be  several 
times  greater  than  those  invacuo,  the  amount  of  effort  would 
still  be  low  enough  for  preliminary  optimization  studies, 
since  the  invacuo  studies  were  not  particularly  computation- 
ally demanding  [12],  particularly  if  the  one  degree  of  free- 
dom equation  yields  satisfactory  accuracy. 
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CHAPTER  5 

RESULTS  AND  COMPARISONS 

Five  example  shells  were  analyzed  in  vacuum  and  in  water 
so  as  to  examine  the  effect  of  immersion  in  water,  examine  the 
accuracy  of  the  formulation  described  herein,  and  compare  the 
accuracy  of  the  Donnell  type  and  Flugge  type  shell  theories. 

In  addition,  for  several  examples  the  effect  of  stiffener  ec- 
centricity is  examined. 

The  first  group  consists  of  three  shells  developed  in  the 
studies  of  Ref  [11]  and  design  examples  using  a minimum  fre- 
quency constraint  from  Ref  [12]  immersed  in  sea  water  at  a 
depth  of  1,000  feet.  They  will  be  referred  to  as  examples  1, 

2 and  3,  respectively.  These  shell  are  analyzed  in  water  and 
invacuo,  where  for  the  latter  the  effects  of  the  hydrostatic 
pressure  of  immersion  are  considered.  Shells  with  outside 
stiffeners  are  also  analyzed  in  addition  to  the  inside  stif- 
fened configurations  of  Refs  [11,12].  The  parameters  used 
for  these  studies  are 

R = 198  in.,  L = 594  in.,  y = 0.33,  E = 30  x 106psi, 

Pc  = 7.33  x 10  4 slug/in3 

Pw=  0.969  x 10  4 slug/in3,  c - 60,000  in/sec,  Depth 

of  immersion  1000  feet 
Where  for  example  1 

h = 1.2056  in.,  ; = 30.17  in.,  h^  = 11.02  in., 

h-2  = 0.3071  in.,  h^  = 10.363  in.,  h^  = 


0 . 2 3 ” 3 in . 
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for  example  2 

h = 1.2216  in.,  ^ = 33.85  in.,  = 20.72  in.,  h,,  = 0.4653  in. 

h_  = 17.55  in.,  h_^  = 0.3950  in. 
and  for  example  3 

h = 1.1623  in.,  i = 16.71  in.,  = 9.411  in., 

= 0.1832  in.,  = 5.937  in.,  = 0.2010  in. 

These  shells  differ  primarily  in  the  nature  of  their  stiffeners. 
Example  2 employes  relatively  large,  widely  spaced  siffeners, 
while  Example  3 uses  relatively  small,  closely  spaced  stiffeners. 

The  effect  of  immersion  in  water,  of  ignoring  the  u and 
v motion  components  and  stiffener  eccentricity  are  very  similar 
for  all  these  examples.  Thus,  only  the  results  of  Example  1 
are  presented  in  tabular  form.  It  may  be  seen  from  Table  1 
and  Figure  5.1  that  the  natural  frequency  characteristics  of 
the  invacuo  model  for  shells  immersed  in  water  is  grossly  in- 
accurate, particularly  with  respect  to  frequency  separation 
which  is  much  less  in  water.  Thus,  the  invacuo  model  does  not 
seem  appropriate  for  optimal  frequency  separation  studies. 

Other  differences  are  also  apparent.  Although  there  is 
a substantial  difference  between  frequencies  computed  using 
three  degrees  of  freedom  [Eq.  (4.22)]  and  one  degree  of  free- 
dom [Eq.  (4.25)]  with  the  invacuo  model,  particularly  where 
n = 1,  this  difference  is  for  practical  purposes  negligible 
for  shells  considering  fluid  interaction  effects  of  water.  The 
largest  difference  for  shells  immersed  in  water  also  occurred 
n = 1.  However,  here  the  difference  was  typically  only  about 


Natural  Frequency  Hz 


Figure,5-I  Natural  Frequencies  of  Example  I Shell 


In  vacuum  In  seawater 


Table  1.  Natural  Frequencies  Hz  Example  1-Inside  Frames 


m 

1 

1 

3 

n 

(1) 

(2) 

mm 

WW 

mm 

WBM 

0 

29.4 

29.4 

39 . 5 

39.5 

46.8 

46.8 

1 

23.0 

22.4 

43.9 

n 

56.7 

56.5 

? 

u 

10.3 

10.2 

29.3 

28.9 

44.6 

44.4 

3 

6.96 

6.90 

19.8 

19.6 

33.4 

33.2 

4 

13.14 

13.08 

18.8 

18 . 7 

28.1 

28.0 

5 

24.8 

24.7 

27.2 

27  . 1 

32.1 

32.0 

6 

40 . 4 

40.3 

■ | 

44.6 

44.5 

0 

162 

158 

159 

158 

158 

157 

1 

107 

76.8 

123 

142 

m 

2 

34.3 

29.5 

81.7 

104 

3 

18.6 

17.6 

49.1 

46.6 

76 . 1 

73.4 

4 

30.8 

30.0 

60.0 

58.6 

5 

53.1 

52.2 

56.8 

55.9 

64.5 

63.6 

6 

80.8 

80.0 

82.2 

85.2 

84.4 

(1)  Three  degrees  of  freedom 

(2)  One  degree  of  freedom 
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3 per  cent.  Thus,  only  one  degree  of  freedom  need  be  con- 
sidered for  such  shells. 

Comparing  Tables  1 and  2,  it  may  be  seen  that  stif- 
fener eccentricity  has  little  effect  except  at  the  lower 
frequencies  where  outside  stiffeners  produced  lower  values. 
Thus,  for  optimization  under  a minimum  frequency  constraint 
inside  stiffeners  seem  more  effective. 

A comparison  of  Donnell  vs.  Flugge  theory  for  Examples 
1 and  2 is  given  in  Table  3.  It  may  be  seen  that  Donnell 
theory  yields  substantially  higher  values  for  the  lower  in- 
vacuo  frequencies.  As  will  be  seen  from  Example  5,  this 
difference  represents  a greater  inaccuracy  on  the  part  of 
the  Donnell  theory  at  these  frequencies. 

Example  4 is  used  to  provide  some  correlation  between 
the  present  analysis  and  experimental  results  for  shells  im- 
mersed in  water.  Unfortunately,  the  differences  in  end  con- 
ditions cloud  the  comparison  somewhat.  Still,  this  com- 
parison is  useful  since  for  such  shells  the  effect  of  end 
condition  is  not  great.  A comparison  with  experiments  in 
air  is  also  provided. 

For  this  example; 

R = 20.25  in.,  L = 60.75  in.,  y = 0.3,  E = 30  x 106psi, 
Pc  = 7.77  x 10 ' 4 s lug/ in3,  pw  = 0.934  x 10  4slug/in3, 
c = 60,000  in/sec,  Depth  of  Immersion  5 feet, 
h = .177  in.,  i = 3.375  in.,  h^  = 1.783  in., 
h7  = 0,  h3  = 0,  h^  = 


0.181  in. 


In  vacuum  In  seawater 


Table  2.  Natural  Frequencies  H:  Example  1-Outside  Frames 


m 

1. 

2 

3.  i 

n 

(1) 

(2) 

mm 

— 

mm 

(2) 

0 

29.4 

29.4 

39.5 

39.5 

46.8 

46.8 

1 

23.2 

22.5 

44  . 3 

43.1 

mm 

57.1 

2 

10.2 

10 . 0 

30.3 

28.8 

46.4 

46.0 

3 

5.85 

5.79 

20.5 

3 5.7 

4 

12.1 

12  . 1 

18 . 4 

in 

mm 

5 

24.2 

24.0 

26.2 

26.0 

32.8 

32.6 

6 

40.0 

39.8 

40.8 

40.6 

44 . 2 

43.9 

0 

162 

158 

159 

158 

158 

157 

1 

109 

77.4 

■MB 

in 

144 

141 

2 

34.5 

29 . 2 

85.3 

77.4 

B9 

108 

3 

15.8 

14.2 

mm 

48 . 2 

82.7 

79.1 

4 

29.0 

27.8 

42.0 

40.2 

65.3 

62.9 

5 

52 . 3 

50.8 

55.2 

53.6 

66.7 

64 . 8 

6 

80.7 

79.0 

80.9 

79.1 

85.3 

83.4 

— 1 

(1)  Three  degrees  of  freedom  analysis 

(2)  One  degree  of  freedom  analysis 


Table  3.  Comparison  of  Frequencies  with  Method  of  Ref  [11] 

Example  1 Example  2 


(1) 

(2) 

(3) 

(1) 

(2) 

(3) 

“1  • 

Hz 

12.03 

17.6 

18.62 

28 . 37 

36.70 

41.50 

(i)  2 ~ 

Hz 

22.90 

29.52 

30.81 

51.96 

58.58 

60.88 

"3  " 

Hz 

30.30 

30.00 

34 .27 

51.96 

68.17 

70.55 

CD 

Bronowicki  et  al 

freedom 

. Flu 

0ge  theory 

three 

degrees 

of 

(2) 

Present 

freedom 

analysis 

Donnell  theory 

one  degree  of 

(3) 

Present 

freedom 

analysis 

Donnell  theory 

three 

degrees 

of 

i 


1 


[ 


— 
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It  may  be  seen  from  Table  4 that  the  one  degree  of  free- 
dom model  used  here  provides  reasonable  agreement  with  experi- 
ment. Considering  three  degrees  of  freedom  produces  slightly 
poorer  results  due  to  the  inaccuracy  inherent  in  the  Donnell 
theory  for  shells  of  such  parameters.  Agreement  is  better 
for  shells  in  water  than  in  air.  This  is,  of  course,  expected 
since  the  inaccuracies  in  the  Donnell  theory  affect  only  the 
shell  stiffener  terms  and  thus  this  improvement  in  accuracy 
in  water  indicates  that  the  fluid-structure  interaction  model 
yields  good  results. 

The  final  example  is  used  to  provide  a more  extended  com- 
parison between  Donnell  and  Flugge  theories.  " The  paremeters 
for  this  study  are: 

R = 10.2  in.,  L = 30.0  in.,  y = 0.35,  E = 10  x 106psi, 

Pc  = 2.54  x 10  4 slug/in3,  pw  = 9.69  x 10 " Sslug/in3, 
c = 60,000  in/sec.,  Depth  of  Immersion  0 feet, 
h = 0.330  in.,  \ = 5.00  in.,  h^  = 1.00  in.,  h?  = h^  = 0, 
h4  = 0.375  in. 

The  results  for  the  invacuo  study  are  compared  in  Table 
5 with  those  given  by  Harari  and  Baron  [4].  It  may  be  seen 
that  at  the  lower  frequencies  the  difference  is  again  quite 
pronounced,  although  proportionally  less  than  in  Examples  1 
and  2,  where  stability  effects  contribute  substantially  and 


thus  Donnell  inaccuracies  are  compounded.  Divergence  between 
theories  increases  with  increasing  m and  decreases  with  in- 
creasing n.  Unfortunately,  at  larger  n,  where  agreement  be- 


In  water 


(1)  Experimental  results  from  [10] 

(2)  Present  one  degree  of  freedom  analysis 


Table  4.  Comparison  with  Paisley  et  al.  [10] 


Natural  Frequencies  Hz 

2 


(1) 

181 

CD 

128/134 

147 

414 

204 

218 

382 

373 

395 

466 

599 

657 

655 

860 

999 

902 

262/270 


369/383 


650 


990 


719 


1127 


760 


1036 


738 


809 


1089 


(1)  Experimental  results  from  [10] 

(2)  Present  one  degree  of  freedom  analysis 
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In  Vacuum 


Table  5.  Comparison  with  Experiment  and  Haran  and  Baron 

Natural  Frequencies  Hz 


CD  (D  (3) 


1232  1133  1526 


640  756 


(2)  (3) 


3020  3105  3 


3 

787 

832 

922 

1190 

1194 

1381 

1602 

1650 

1923 

m 

1310 

1432 

1483 

1503 

1575 

1764 

1806 

1826 

2183 

5 

1938 

2253 

2253 

2059 

2331 

2483 

2276 

2474 

2839 

6 2594  3276  1 3205 


7 3179  4466 


8 3728  5853  5637 


2602  3424  3755 


(1)  Experimental  results  Ref  [18] 

(2)  Harari  and  Baron  orthotropic  theory  [4] 

(3)  Present  one  degree  of  freedom  analysis 
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tween  theories  is  good,  orthotropic  theory  produces  poor  re- 
sults [ 4 ] . 

Thus  it  appears  that,  although  the  present  method  gives 
reasonably  good  results  for  shells  immersed  in  water,  use  of 
the  Flugge  theory  could  significantly  improve  accuracy. 

Example  4 was  also  examined  for  effects  of  stiffener 
eccentricity  and  ignoring  the  u and  v motions  as  was  Example 
5,  which  was  also  analyzed  for  frequencies  immersed  in  water. 
Results  produced  followed  closely  those  of  Examples  1-3  with 
regard  to  these  effects  except  that  the  effect  of  stiffener 
placement  and  use  of  only  one  degree  of  freedom  were  even  less 
significant  in  Examples  4 and  5. 


CHAPTER  6 


CONCLUSION 

The  interim  procedure  for  calculation  of  frequencies  of 
submerged  shells  described  here  seems  quite  adequate  for  pre- 
liminary optimization  studies  of  shells  immersed  in  water. 

The  simplicity  of  the  shell  model  and  the  fact  that  only  one 
degree  of  freedom  need  be  considered  allows  rapid  evaluation 
of  frequency  behavior,  an  important  asset  for  optimization, 
which  usually  requires  several  hundred  sets  of  such  evalua- 
tions for  the  determination  of  one  optimal  design.  For  all 
modes  of  interest  here  the  radiated  pressure  was  real  and  thus 
the  procedure  readily  allows  consideration  of  the  fluid  inter- 
action effect.  Thus,  the  existing  method  seems  well  adapted 
for  preliminary  optimization  under  a minimum  frequency  con- 
straint or  for  preliminary  optimal  frequency  separation  prob- 
lems involving  the  lower  frequencies. 

Use  of  a Flugge-type  shell  theory  would  improve  accuracy 
and  make  the  procedure  useful  for  design  purposes  for  such 
problems.  Extension  to  treat  cases  where  the  radiated  pressure 
is  complex  is  also  straightforward. 

The  primary  limitation  of  this  method  even  with  the 
above  extensions  in  treating  the  optimal  frequency  separation 
problem  is  the  possible  presence  of  active  inter-ring  vibra- 
tion modes  at  the  optimum  [13].  Realistic  treatment  in  such 
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circumstances  requires,  use  of  discrete  theory  [4],  [19]  and  its 
associated  large  computational  effort.  Even  this  extension  is 
probably  tractable  in  light  of  current  optimization  develop- 
ments [20,21].  The  other  important  limitation  results  from 
the  fact  that  the  end  effects  of  fluid  interaction  with  the 
shell  are  ignored.  These  effects  would  play  an  important  role 
in  frequencies  associated  with  primarily  axial  models,  although 
these  fluid  effects  should  be  less  important  here  than  they  are 
for  primarily  axial  modes.  The  inclusion  of  these  effects  is 
a formidable  problem  particularly  if  an  approach  considering 
these  effects  is  to  be  used  for  the  purposes  of  optimization. 
Thus,  the  prospect  of  early  development  of  reasonably  accurate 
optimization  procedures  for  primarily  axial  frequency  separation 
seems  remote. 
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